Electronic Detection of Collective Modes of an Ultracold Plasma 
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Using a new technique to directly detect current induced on a nearby electrode, we measure plasma 
oscillations in ultracold plasmas, which are influenced by the inhomogeneous and time- varying den- 
sity and changing neutrality. Electronic detection avoids heating and evaporation dynamics asso- 
ciated with previous measurements and allows us to test the importance of the plasma neutrality. 
We apply dc and pulsed electric fields to control the electron loss rate and find that the charge 
imbalance of the plasma has a significant effect on the resonant frequency, in excellent agreement 
with recent predictions suggesting coupling to an edge mode. 
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Collective oscillations are central to the study of plas- 
mas as they embody the rich physics unique to the 
plasma state and provide diagnostics of plasma density 
and temperature. Ultracold neutral plasmas (UCPs) 
PQ [2] are a novel system for the study of collective behav- 
ior as they have extremely low temperatures (1-100 K), 
an inhomogeneous density and are on the border of strong 
coupling. Previous experiments on UCPs have observed 
plasma oscillations [3 j, Tonks-Dattner resonances [4], ion 
acoustic waves [5j [6] and a high-frequency electron drift 
instability [7]. 

Like most laser-produced plasmas, UCPs are freely ex- 
panding systems. External fields can force varying rates 
of electron evaporation, which gives a time-dependence to 
the global neutrality, an important feature of these sys- 
tems that has only recently been explored [8] [9] . A recent 
theoretical study of plasma oscillations [9 , the simplest 
collective mode, considered the impact of electron loss 
and predicted the existence of a zero-temperature mode 
with a resonant frequency that increases for less neutral 
plasmas. The theory assumed a spherically symmetric 
electron spatial distribution, n e , but we have found a 
significant asymmetry in typical experiments [8J and the 
effect of this on electron oscillations has not been ad- 
dressed. 

Here, we excite and detect plasma oscillations of an 
UCP taking into account the changing density, neutral- 
ity and symmetry of n e . We find the variation of the 
resonant frequency with neutrality agrees well with the 
predictions of [9]. In addition, we have developed a new 
diagnostic for UCPs where we directly detect oscillations 
through the current induced on a nearby electrode. The 
method is more accurate than previous measurements 
based on electron emission, can be used when charged 
particle detection is not possible, and gives UCP experi- 
ments access to a broader spectrum of collective modes. 

We create a plasma by two-photon ionization of about 
10 6 metastable Xe atoms collected in a magneto-optical 
trap [10] . The initial plasma density is roughly Gaussian 
with an rms radius of 0.3 - 0.6 mm. The initial energy 



given to the electrons, E e , is controlled by tuning the 
energy of the ionization laser above the ionization limit. 
After creation, the plasma loses a few percent of the elec- 
trons until a sufficient charge imbalance exists to trap the 
remaining electrons, forming a plasma. The plasma ex- 
pands, driven by thermal electron pressure. The ion dis- 
tribution, rii, remains Gaussian, following a self-similar 
expansion [11] described by of (t) = of (0) + v^t 2 , where 
<Ji is the rms radius of the ion distribution and the expan- 
sion velocity is typically = 50 — 100 m/s, determined 
by£ e . 

Two wire mesh grids located 1.4 cm on either side of 
the plasma apply a weak electric field (5-10 mV/cm) that 
directs evaporating electrons out of the plasma region 
and onto a microchannel plate detector. A typical elec- 
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FIG. 1. Comparison of rf signals to electron signals, (a) The 
dotted and solid black lines are the electron emission signals 
with and without a / = o; r f/27r = 14 MHz rf field applied. The 
rf phase and amplitude signals are the difference of the signal 
obtained with and without the plasma, have been rescaled by 
the same arbitrary factor and are offset for clarity, (b)-(d) 
Rf phase change compared to subtracted electron signal for 
three different frequencies. 
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tron current signal including the prompt loss of electrons 
and the electron evaporation during expansion is shown 
in Fig.^a). 

UCPs are small systems, consisting of only 10 4 — 10 9 
ions and electrons and sizes of 0.1 mm to 1 cm, so avail- 
able experimental probes have been limited. Optical 
absorption and fluorescence imaging of the plasma ions 
[12] [13] have provided spatially resolved density and ve- 
locity measurements of the ions. Information about the 
electrons has predominantly been obtained by monitor- 
ing the loss rate on a charged particle detector as de- 
scribed. This method of electron detection has succeeded 
in observing plasma oscillations [3] [4] , but only indirectly 
by applying a constant driving field that resonantly heats 
the electrons and observing an enhanced loss rate. Thus 
the measurements are subject to the dynamics involved 
with heating the electrons and their subsequent evapora- 
tion. 

We present a new approach to studying electron reso- 
nances in UCPs by directly measuring changes in the rf 
field. Measurements of rf absorption are commonplace 
in low density laboratory plasmas. Most analogous to 
our system are measurements of zero-temperature oscil- 
lations in non-neutral plasmas trapped in Penning traps 
[T4HT7] . These plasmas are typically of similar size and 
density to our neutral plasmas, but our rapidly expand- 
ing plasmas are untrapped so the measurements must be 
made during the fast time evolution of the plasma den- 
sity. Our resonances last only a few /is, about 100 times 
shorter than the averaging times used for non-neutral 
plasmas. 

We detect plasma modes by applying a weak, contin- 
uous rf drive at frequency cj r f to the grid located above 
the plasma and monitor the amplitude and phase changes 
of the voltage coupled to the symmetric grid below, V5, 
as sketched in Fig. [2] In the absence of a plasma, Vb 
is constant in time and simply related to the electrode 
geometry. When a plasma is present and driven near a 
resonant frequency, the oscillation of the plasma induces 
a current on the bottom grid that interferes with the 
background signal. The signal from the plasma is much 
smaller than the background, so we observe only small 
changes in the amplitude and phase of V& as the plasma 
density quickly scans through resonance with the driving 
field. The results of both quadratures of this homodyne 
measurement are shown in Fig. [TJa). All measurements 
are done with > 400 kHz bandwidth, sufficiently large to 
capture the fast changes in the rf signals. Due to this 
large bandwidth and small particle number, the signal- 
to-noise ratio on a single experimental shot is often less 
than 1, so all plots are an average of at least 150 shots. 

We qualitatively reproduce the shapes of the rf signals 
using a circuit model that treats the plasma as a series 
RLC oscillator in close analogy with [H] . This shows 
that a peak in the phase-change signal corresponds to 
the time at which the plasma is resonant with the ap- 
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FIG. 2. Schematic for homodyne detection of plasma oscilla- 
tions. The rf source is attenuated (Atten.) and applied to the 
top grid. The dc bias and voltage pulses are applied at Vb. 



plied field, and thus we will focus on this signal in the 
remainder of this paper. A more detailed analysis could 
also provide information about mode damping and the 
electron temperature. 

Figure [l] shows the comparison of the rf measurements 
with the previously employed enhanced electron emission 
method. The earliest-time peak in both the rf and elec- 
tron signals is the zero-temperature plasma resonance. 
For large plasmas and higher cj r f, we routinely see ex- 
tra peaks at later time in the electron signal, e.g. three 
peaks in parts (a) and (b) of Fig. [I] that were previously 
identified as temperature-dependent Tonks-Dattner res- 
onances [4]. We have difficulty resolving more than one 
extra peak in the rf signals. For larger driving amplitude 
we can resolve a total of three peaks in the rf signals, but 
the third is a factor of 10 or more smaller in amplitude 
than the first. We anticipate that a higher signal-to-noise 
measurement would allow us to see more of the modes. 
Molecular dynamics simulations [18] have also found mul- 
tiple peaks in the rf absorption in UCPs, but a detailed 
understanding of these modes is still lacking. We focus 
our analysis on the zero-temperature resonance. 

The resonance times in the two methods agree to 
within 1 /is in most cases, after correcting for a delay 
due to finite bandwidth. There is disagreement at low u r f 
(< 6 MHz), where the electron signal becomes less reli- 
able. The electron signal measures evaporation caused by 
rf heating, but we should not expect a linear relationship 
between energy absorption and electron emission. As the 
heating begins and electrons are lost, the plasma poten- 
tial well deepens, so a greater input energy is needed 
to subsequently maintain the same electron flux. This 
effect is most evident at low frequencies as the plasma 
response becomes much broader in time. By contrast, 
the rf signals measure an induced current that is directly 
proportional to the amplitude of electron oscillations. 

To understand the cold plasma resonance, we consider 
the full picture of the spatial distribution of electrons. 
Optical measurements of the plasma ions have shown 
that rii remains Gaussian throughout expansion [TT]. At 
the center of the plasma, n e must be nearly equal to n^, 
but electron loss ensures deviations at the plasma edge. 
Bergeson and Spencer [19] solved the cold plasma fluid 
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FIG. 3. Edge- mode theory compared to free expansion and 
pulsed emission data. The edge-mode and quasi-mode rem- 
nant curves are from [9]. Blue points with error bars are the 
average of many free expansion measurements. Red points 
are individual measurements after pulsed electron emission. 
Results are the combination of many experimental runs all 
with E e /k b = 100 K and 6 x 10 5 < N t < 10 6 . 



FIG. 4. Electron distributions with 5 = 0.53. The shaded re- 
gion represents the Gaussian ion distribution. Cold electron 
distributions match m inside the boundary lines and drop 
sharply to zero outside. The solid line is the symmetric elec- 
tron boundary from theory. The dashed line is the electron 
boundary during free expansion and the dotted line is the 
boundary after a pulsed emission of electrons, both of which 
have an asymmetry caused by the external field, E ex t- 



equations for a perfectly neutral plasma (n e = Ui in all 
space), assuming no electron loss, and found only a sin- 
gle quasi-mode with a maximum energy absorption at 
a frequency uo = 0.24co> p o, where uj p q = ^/e 2 n e o/m e eo 
and n e o is the central plasma density. But even early in 
the plasma lifetime the deviation from neutrality may be 
non- negligible, owing to the prompt loss of electrons at 
plasma creation [Fig. Qa)]. Lyubonko et al. [9] allowed 
for electron loss by treating ?7/g as a truncated Gaussian 
with the truncation radius set by the charge imbalance 
5 = (Ni — N e )/Ni, where Ni and N e are the number 
of ions and electrons. For plasmas with any significant 
charge imbalance (5 > 5%), they found a mode where the 
majority of energy was absorbed near the sharp edge of 
the cold electron distribution. The relative frequency of 
the edge- mode resonance, uo re i = uj/uj p o, increases with S 
as shown in Fig. [3J 

We perform two types of experiments to test the edge- 
mode prediction. First, we record the resonance times 
from the rf phase change signal during normal free expan- 
sion of the plasma. Figure [I] shows typical measurements 
for different cj r f . For smaller cj r f , the resonances are ob- 
served later in time, corresponding to a lower density. 
But due to the continuous electron loss, later time also 
corresponds to larger charge imbalance, which should in- 
crease the relative resonant frequency. We fit a Gaussian 
to the first peak of the rf phase change to get the reso- 
nance time. We independently measure the plasma ex- 
pansion velocity and ion number, which allows us to cal- 
culate the relative resonant frequency uj re \. The charge 
imbalance is calculated by integrating the electron sig- 
nals. Only a small 10-15 mV/cm electric field is needed 
to collect all plasma electrons on our detector. The frac- 
tion of the integrated current that arrives before the res- 
onance time gives the charge imbalance S. The results 



are shown in Fig. [3] 

A potential problem with this test is that we expect 
n e in the experiment to have a different shape than the 
symmetric truncated Gaussian used in theory. We have 
found that n e in a freely expanding plasma can develop 
a strong asymmetry from the influence of externally ap- 
plied or stray dc electric fields, and we intentionally ap- 
ply such a field to facilitate electron detection. This field 
perturbs n e , which we have used to explain the observed 
rate of electron loss from our system [8]. The difference 
between our calculated n e and the symmetric case for the 
same charge imbalance is shown in Fig. [4j 

In a second experiment, we dump electrons from the 
plasma with short voltage pulses of 0.5-2 /xs, chosen to be 
longer than the electron collision time. This creates an 
n e (dotted line in Fig. [4| that is closer to the symmetric 
distribution used in theory and gives more control over 
5. Fig. [5] shows examples of the electron emission signal 
with pulses applied. 

Before the voltage pulse, we assume a Gaussian and 
an asymmetric n e , like the dashed line in Fig. |4j After 
the pulse, the remaining electrons will be concentrated 
mostly at the center of n^. The dc electric field will still 
polarize the plasma, but the electrons are now held in a 
deep well and there are not enough to reach the edge of 
the ions. Electron emission ceases for 10 to 20 [is. As the 
plasma continues to expand, the well depth decreases and 
n e slowly becomes less symmetric until it again reaches 
the edge of ni and electron emission returns. We can 
calculate n e after a pulse using the same algorithm used 
in [8] but fixing the value of 5. Fig. [4] shows one example 
result. 

The drop in N e must affect the ion expansion, as the 
ions outside of the electron cloud are free to move in 
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FIG. 5. Electron signals (left) and rf phase change signals 
(right) with a voltage pulse of 1 fis length and varying ampli- 
tude applied at 40 /xs. The applied frequency is / = 4 MHz. 



response to both the external field and their own inter- 
actions without electron screening. For larger electron 
dumps, the density of free ions is large enough that a 
significant Coulomb explosion occurs and disrupts the 
following neutral plasma expansion. The electron signal 
following a pulsed emission is an indicator of the effect 
of this Coulomb explosion on the plasma. Since the elec- 
tron loss rate is related to the size and shape of the ion 
cloud, an electron signal very similar to that during nor- 
mal free expansion indicates a minimal change in plasma 
expansion. As shown in Fig. [5j we can for lower voltage 
pulses get excellent agreement in the late-time electron 
signals. 

For a given voltage pulse, we adjust co> r f and look for 
resonances in the rf phase change during the dead time of 
electron emission, as shown in Fig. [5] When we are near 
a resonance, a clear peak is observed that we identify 
as the edge-mode resonance. Changing the strength of 
the voltage pulse dumps more electrons, increasing <$, but 
leaves all other parameters unchanged. We see that this 
direct change in S increases the time of resonance, which 
corresponds to an increase in u re i as predicted. 

The data in Fig. [3] are created by applying voltage 
pulses at many times during expansion. At each time, we 
sweep o; r f to find resonance peaks that come just after the 
pulsed emission. We calculate uj re \ from the resonance 
time assuming an uninterrupted ion expansion, and we 
only plot points where the observed peak comes within 
7 /is of the end of the voltage pulse. To get to higher 
S we wait until later time in the expansion, letting the 
plasma naturally lose more charge, before applying the 
pulse. During pulsed emission, we pulse the electrons 
away from the detector but still collect all other electrons. 
We find S by comparing the integrated current after the 
time of resonance to the full integrated current in the free 
expansion experiment. 

It is clear from Fig. [3] that the free expansion and 
pulsed experiments agree with each other and the edge- 
mode theory. Agreement in the free expansion data is 
worst at large (5, but these rf signals are very broad, span- 
ning 10s of /is. Significant changes to collision rates and 
significant electron loss during the response time may af- 



fect the magnitude of the rf signal and distort the results. 
The pulsed emission experiments give us control over 5, 
and we can more easily probe large values. It is unclear 
why the rf signal after a pulsed emission has a stronger 
and sharper response than during free expansion. 

The importance of the shape of n e is also not imme- 
diately clear. During free expansion, we expect n e to 
be more asymmetric than the distribution after a pulsed 
emission, yet both measurements seem to agree equally 
well with the perfectly symmetric theory, suggesting that 
the spatial distribution of charge is only of secondary im- 
portance to the integrated amount of charge. It seems 
surprising that the position of the electron cloud edge is 
not a larger factor given the finding in [9 that the large 
majority of energy is absorbed at this edge. It would 
be interesting to see how the solution of the cold plasma 
fluid equations changes for the asymmetric distributions 
of Fig. |4j 

In conclusion, we have shown detailed measurements of 
cold plasma oscillations in an expanding ultracold plasma 
with a time- varying n e . Both of our experimental ap- 
proaches support an increase of resonance frequency with 
charge imbalance, in agreement with a zero-temperature 
theory. We have also presented a new diagnostic tool 
for probing oscillations in ultracold plasmas that is more 
accurate and more versatile than the previous method. 

An important feature of our new measurement is that 
it allows us to probe electron properties through their 
resonant behavior without the need for charged particle 
detection. This advantage was evident in our ability to 
observe resonances during the dead time in electron emis- 
sion. Charged particle detection is also prevented when 
a magnetic field is applied transverse to the axis of the 
electric grids. We have observed upper hybrid resonances 
in this setup, which will be the subject of future work. 

We have, in some regimes, also been able to observe the 
free decay of electron oscillations after abruptly turning 
off the rf drive. Studying mode damping should pro- 
vide information on the collision properties and electron 
temperature in UCPs. It is worth noting that we ob- 
serve plasma resonances at much later times than are 
normally studied in UCPs. At 90 /is, the plasma density 
has dropped to 2 x 10 5 cm -3 , but we can still observe 
clear collective behavior, which is a testament to the ex- 
tremely low electron temperature, expected to be less 
than 1 K at that time. 
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